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Implementation of a pressure gradient method of impedance eduction in two NASA Langley flow ducts 
is described. The Grazing Flow Impedance Tube only supports plane-wave sources, while the Curved Duct 
Test Rig supports sources that contain higher-order modes. Multiple exercises are used to validate this new 
impedance eduction method. First, synthesized data for a hard wall insert and a conventional liner mounted 
in the Grazing Flow Impedance Tube are used as input to the two impedance eduction methods, the pressure 
gradient method and a previously validated wall pressure method. Comparisons between the two results are 
excellent. Next, data measured in the Grazing Flow Impedance Tube are used as input to both methods. 
Results from the two methods compare quite favorably for sufficiently low Mach numbers but this comparison 
degrades at Mach 0.5, especially when the hard wall insert is used. Finally, data measured with a hard wall 
insert mounted in the Curved Duct Test Rig are used as input to the pressure gradient method. Significant 
deviation from the known solution is observed, which is believed to be largely due to 3-D effects in this flow 
duct. Potential solutions to this issue are currently being explored. 
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= mode coefficients upstream and downstream of liner test section (Pa) 

= speed of sound (m/s), uniform flow Mach number, and mean static density (kg/m 3 ) 

= transcendental function, source frequency (Hz), and time (s) 

= height, width, and length of duct (m) 

= y/—l, unit imaginary number 
= freespace and axial wavenumbers (m 1 ) 

= axial locations of the leading and trailing edges of liner (m) 

= number of horizontal and vertical hard wall duct modes 
= number of grid points in the x and z directions of grid 
= 3D and quasi-3D acoustic pressure fields (Pa) 

= source and termination plane acoustic pressure profiles (Pa) 

= axial, vertical, and horizontal coordinates (m) 

= mode shape, complex coefficient, and mode eigenvalue (m 1 ) for vertical mode 
= normalized acoustic admittance of lower and upper wall liner 
= normalized acoustic impedance of lower and upper wall liner 
= unconstrained design variables for optimization algorithm 
= normalized acoustic resistance of lower and upper wall liner 
= normalized acoustic reactance of lower and upper wall liner 
= normalized acoustic conductance and susceptance of lower wall liner 
= objective function (Pa/m), circular harmonic frequency (s -1 ), and complex absolute value 
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Superscripts: 

+,— = right and left moving wave indicators 

Subscripts: 

m,n = horizontal and vertical mode orders 

/, / = axial and vertical grid line counters 

I. Introduction 

As part of NASA’s Fundamental Aeronautics Program, focused efforts are currently underway to develop concepts 
to confine subsonic aircraft objectionable noise within the boundaries of the airport. Although acoustic liners mounted 
in the aircraft engine nacelles provide a significant portion of the current fan-noise reduction, they must be further 
optimized if this goal is to be achieved. This will require more detailed understanding of an intrinsic parameter of the 
nacelle liner referred to as the locally-reacting acoustic impedance. This intrinsic parameter is dependent on sound 
pressure level and grazing flow velocity. Even after four decades of research in impedance measurement technology, 
accurate and reliable tools for a direct measurement of this intrinsic parameter have proved challenging. 

As an alternative to direct measurement of the locally-reacting acoustic impedance, various numerical methods 
(which are used in conjunction with measurements from experimental facilities) have been developed. The numerical 
methods 1-12 include codes based on single mode, multi-mode, and finite element analyses. Of these, finite element 
codes developed at the NASA Langley Research Center (LaRC) have been thoroughly tested in the open duct acoustics 
literature 2 * ' 7 ' 9 and have been demonstrated to provide quality results when used in conjunction with measurements 
acquired in the NASA LaRC Grazing Flow Impedance Tube (GFIT). A schematic of the GFIT’s flow duct is given in 
Fig. 1. The figure shows the liner test section (i.e., the portion of the flow duct that contains the liner) along with the 
hard wall sections of duct upstream and downstream of the liner test section. The mean flow is directed from left to 
right along the axis of the duct. The locally-reacting acoustic impedance is educed for test liners in the GFIT using a 
procedure that has been well established. This procedure 2 uses the acoustic pressure measured by microphones flush- 
mounted in the hard wall opposite the test liner (see Fig. 1) to educe the test liner impedance and has been referred to 
in subsequent works as the wall pressure method (WPM) of impedance eduction. 

A recent paper 13 gives a description of the design and development of the Curved Duct Test Rig (CDTR). This rig 
was initially designed as a test bed to quantify the performance of current liner designs and to develop new concepts 
for efficient duct liner treatments in the presence of grazing flow and wall curvature. Although impedance eduction 
at LaRC is currently performed using acoustic measurements acquired in the GFIT, the decision was made in 2010 to 
acquire the capability to perform impedance eduction using measurements obtained in the CDTR. This decision was 
made because: (1) flow profiles acquired in the CDTR are typically much closer to those of a full-scale engine than 
those acquired in the much smaller GFIT, (2) the CDTR waveguide is closer to the size of the aft bypass duct of a 
business jet, and (3) higher-order modes can be controlled in the CDTR. Figure 2 gives a schematic of the flow duct for 
the CDTR with straight walls. Shown in the figure are the liner test section (rotated ninety degrees for convenience) 
and the hard wall sections upstream and downstream of the liner test section. Note that the liner test section contains 
test liners on opposite walls and is not populated with flush mounted microphones as in the GFIT. 

In a previous paper, 14 it was demonstrated that the GFIT impedance eduction method could be successful if the 
impedance eduction data were acquired in a hard wall adjacent to the test liner (i.e., the left or right side wall in Fig. 1). 
While plans are to investigate this method of impedance eduction in the CDTR, the hardware required for obtaining the 
necessary sidewall measurements is not currently in place. However, in a recent paper, 15 a new numerical procedure 
for educing the impedance of test liners located in a duct with grazing flow was presented. This procedure, referred to 
as the pressure gradient method (PGM) of impedance eduction has the advantage that it is applicable to ducts with test 
liners on more than one wall as illustrated in Fig. 2, is less cumbersome than the WPM (i.e., no microphones need be 
installed in the liner test section), and can be applied using the instrumentation already installed in the CDTR. Further, 
in Ref. 15 the PGM was successfully tested in a duct with small cross-sectional dimensions (i.e., the GFIT) using 
synthesized data for which only plane waves are cut on in the hard wall sections of the duct. Whether the PGM can be 
successfully applied in a duct with higher-order modes cut on or using measured data remains an open question. 

The current paper has four primary purposes: (1) to describe the implementation and validation of the PGM to 
a duct with small cross-sectional dimensions (i.e., the GFIT), (2) to describe the implementation of the PGM to a 
much larger duct (i.e., the CDTR) for which higher-order modes are cut on, (3) to present a consistency check on the 
implementation of the PGM in the larger duct by demonstrating that it reproduces the same impedance as measured in 
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the GFIT when synthesized data is applied, and (4) to present the results of a validation exercise that uses measured 
data. The remainder of this paper is organized as follows. Section II presents the governing differential equations and 
boundary conditions. Section III provides a description of the pressure gradient method (PGM) of impedance eduction. 
Section IV describes how the PGM is implemented in the GFIT and CDTR, results are presented in section V, and the 
primary conclusions relevant to the study are provided in section VI. 


II. Governing Equations and Boundary Conditions 


Schematics of the flow duct for the GFIT and the CDTR (rotated ninety degrees for convenience) are given in 
Figs. 1 and 2, respectively. Here, the mean flow is directed from left to right along the axis of the duct and for purposes 
of this investigation there is no wall curvature in the CDTR (Fig. 2). The same orthogonal coordinate system is used to 
describe both flow ducts in which z is the axial coordinate, y is the horizontal coordinate, x is the vertical coordinate, 
and the source and termination planes (of each flow duct) are located at z = 0 and z = L, respectively. In the current 
setup, the test liner for the GFIT is located on the upper wall, whereas the CDTR may contain test liners on the upper 
and lower walls. Each liner is assumed to be locally reacting with uniform impedance. The left and right sidewalls 
are rigid in both duct configurations with the leading and trailing edges of each liner located at z = L\ and z = La, 
respectively. 

The analysis used in this paper makes three simplifying assumptions: (1) the flow field in the duct can be decom- 
posed into an acoustic and mean flow field and the mean velocity field has only a uniform axial component, (2) the 
acoustic field is linear and time-periodic, and (3) the sidewalls of the flow duct are rigid. Under these simplifying 
assumptions, the acoustic pressure field in the flow duct can be decomposed into a series of hard wall duct modes of 
the form 

P(z,x,y) = £ cos (“Er) P™( z ’ x ) (!) 

m = 0 

where a time dependence of the form e mt has been assumed. The quasi-3D acoustic pressure field, p m {z,x ), in Eq. (1) 
can be shown to satisfy the convected Helmholtz equation 
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Note that Eq. (2) neglects the effects of the mean boundary layer and will become less accurate at the higher flow 
speeds where gradient in the mean boundary layer are more important. Subsequently, the accuracy of the impedances 
educed using the method proposed in this paper is expected to degrade at the higher flow speed. 

The quasi-3D convected Helmholtz equation (Eq. (2)) is solved subject to the sound source 


Pm{0,x) = S m (x) 


(3) 


and termination plane boundary condition 


Pm(L,x) = Q,„(x) 


(4) 


Finally, each acoustically treated wall of the duct will require an impedance boundary condition. The boundary 
condition for a locally-reacting impedance wall has been derived in a previous paper. 16 When written along the top 
wall the locally-reacting wall impedance boundary condition is 
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and along the lower wall the boundary condition is 
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In Eqs. (5) and (6), i^o and Cn are the unknown normalized impedances of the test liners located on the lower and 
upper wall, respectively. Further, the reciprocal of the normalized impedances are taken as zero along the hard wall 
portion of the duct wall upstream and downstream of the liner test section. Throughout this paper all impedances are 
normalized by dividing by the characteristic impedance, poco, of the air in the duct. 
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For specified wall impedances Co and C//, Eqs. (2)-(6) constitute a well-posed boundary value problem (BVP) that 

and 


can be solved to determine a unique set of acoustic pressure gradient profiles ( ^Pm(0,x) 


dz 


-), normal to 


the source and duct termination planes. The uniqueness of these acoustic pressure gradient profiles provides the basis 
for the PGM impedance eduction method that is the subject of the following section. Closed-form solutions for these 
profiles are not yet available for arbitrary source and termination plane acoustic pressure profiles and wall impedances. 
In this paper, a finite element method that employs cubic Hermite polynomial basis functions 2 is used to numerically 
extract these profiles. 


III. The Impedance Eduction Algorithm 

The PGM has been described in an earlier paper. 15 Therefore only enough is presented in this section to provide 
continuity with the previous paper and to further motivate the methods and goals of this investigation. Because the 
PGM uses an optimization algorithm as a workhorse to extract the impedance, the optimization algorithm is briefly 
described. 


A. Brief Description of the PGM 

The PGM may be summarized as follows: 

1. Measure the mean flow field, the source and termination plane acoustic pressure profiles, and the acoustic 
pressure gradient profiles normal to the source and duct termination planes. 

2. Using initial guesses for the test liner impedances, Co and C//, solve the BVP defined by Eqs. (2)-(6) to obtain 
a unique set of acoustic pressure gradient profiles normal to the source and duct termination planes. Because 
exact solutions to the BVP are generally not available, a numerical method is required to obtain these profiles. 

3. Compare the measured acoustic pressure gradient profiles at the source and duct termination planes to those 
obtained from the numerical solution to the BVP. If the numerically computed and measured acoustic pressure 
gradient profiles match (to within some tolerance) then the initial guesses for the test liner impedances are the 
impedances of the test liners. 

4. If the numerically computed and measured acoustic pressure gradient profiles do not match, then the unknown 
impedances, Co and C//, are iteratively updated using an automated gradient-based optimizer until the measured 
and numerically computed acoustic pressure gradient profiles match (to within some tolerance). The resulting 
impedances are the unknown impedances of the test liners. 

It is noted that the automated gradient-based optimizer iterates on the impedances by minimizing the objective 
function, <j) (^o 7 C//), that is defined as 


<KCo,U) 



dp m (0,xj) 

dz 
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dz 


dp m (L,xj) _ dp m (L,xj ) 
dz dz 
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In Eq. (7), ) j s the normal acoustic pressure gradient obtained from the numerical solution to the BVP and is 

the normal acoustic pressure gradient obtained from the measurement. The automated gradient-based optimizer that 
finds the impedances that minimizes the above objective function is one of the most important features of the PGM. 
Therefore, the strengths and weaknesses of the automated optimizer that is used with the PGM are briefly discussed 
in the following subsection. 


B. The Optimization Algorithm 

The optimization algorithm that is used in this paper is a modification of the Davidon-Fletcher-Powell (DFP) opti- 
mization algorithm. The DFP optimization algorithm is a technique for finding the minimum of an objective function 
such as that defined in Eq. (7). It was first proposed by Davidon 17 and was later reformulated and popularized by 
Fletcher and Powell. 18 The version of the algorithm that is used in this paper has been modified to accept numerical 
differentiation of the gradient as suggested by Stewart. 19 This modified algorithm has become quite popular because 
of its speed and robustness and is referred to here as “Stewart’s adaptation of the Davidon-Fletcher-Powell” (SDFP) 
optimization algorithm. 
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The most important assumption made by the SDFP optimization algorithm (that has bearing on results presented 
in this paper) is the assumption that the objective function can be locally approximated as a quadratic in the region 
around the minimum point. This optimization method, therefore, has the disadvantage that it may converge to a local 
minimum and may become stuck in the portion of the impedance plane where the objective function is extremely flat. 
However, this shortcoming of the SDFP algorithm is more than mitigated by the fact that it converges faster and yields 
more accurate results than many of its competitors. 

A second important issue that affects the accuracy of the SDFP optimization algorithm is that it accepts only real 
design variables and is designed for unconstrained optimization problems. The normalized impedance of each liner is 
therefore divided into its normalized resistance (0o, 0//) and normalized reactance components (%o ,%h) 

Co = 9o + 0 Co (8) 

Cff = 6// + ilH (9) 

Because the optimization algorithm is unconstrained, the SDFP optimization algorithm allows the design variables to 
take on values corresponding to any real number, — °° < (0o, 0//) < < (%o ,%h) < °°- From a practical consid- 

eration of the geometry of nacelle liners used in aircraft engine noise applications (i.e., hole size, face sheet thickness, 
open area ratio, depth of the honeycomb core, etc.), the normalized resistance and reactance of the liner fall within the 
range, 0 < (0o, 0//) < 10, — 10 < (Xo,3C h) < 10. The capability to constrain both the normalized resistance and reac- 
tance to physically realizable values is implemented in this study. This is done in the following manner. Let the design 
variables used by the SDFP optimization algorithm be denoted by q and q. Now introduce the following constraints 
into SDFP: 

0 = 5 [1 +cos (£)] (10) 

X = 10cos(q) (11) 

Note that as the design variables q and q take on any real number, the normalized resistance and reactance of the liner 
will take on only physically realizable values (i.e., 0 < 0 < 10, — 1 0 < % < 10). Finally, a stopping criterion (i.e., 
tolerance) of 1 x 1 0 s is used to terminate the SDFP optimization algorithm. 


IV. Implementation of the Impedance Eduction Method 


The PGM has been implemented in both the GFIT and the CDTR. Because only plane waves are cut on in the hard 
wall sections of the GFIT while higher order modes are cut on in the CDTR, the implementation of the PGM method is 
different in both rigs. Four acoustic measurements are needed to implement the PGM. The first two measurements are 
the source and termination plane acoustic pressure profiles (S m (x) and Q m (x)). The next two acoustic measurements 


are the source and termination plane acoustic pressure gradient profiles ( 


dp»,(0>*) and ggm(M) 


dz 


dz 


-). The following two 


subsections describe how the PGM is implemented in both rigs to acquire these two sets of data. 


A. Implementation in the GFIT 

In the GFIT, two microphones (not shown in Fig. 1) are flush-mounted in the hard wall section upstream of the liner test 
section to provide a measurement of the acoustic pressure gradient profile at the source plane. These two microphones 
are flush-mounted in the lower wall (i.e., at x = x\ = 0) at points in the finite element grid. The first microphone is 
located at the source plane (i.e., z = Z\ =0) and the second microphone is located one grid point to the right of the 
source plane (i.e., z = Z2 )• Therefore the measured acoustic pressure gradient profile normal to the source plane is 
approximated using the forward difference expression 

\Pm(Z2,x) - Pm(zi,x)\ 


dz 


( 12 ) 


(Z2-Zl) 

Because only plane waves are cut on in the GFIT, the acoustic pressure profiles (p m (z 2 ,x) and p m (z\ ,x)) are in- 
dependent of the vertical coordinate, x. Therefore, p m (z 2 ,x) and p m (zi,x) are obtained from the acoustic pressure 
measurements obtained from the two lower wall microphones (i.e., p m (z. 2 ,x) = p m (z 2 ,xi) and p m (zi,x) = p m (z\ ,X| )) 
and only the zeroth order horizontal mode (i.e., m = 0) need be analyzed. 

A method similar to that discussed in the preceding paragraph is used to approximate the acoustic pressure gradient 
profile normal to the duct termination plane. Thus, the backward difference expression 

dgmCM) = \Pm(z nzt x)~ Pm (Znz- 1A')] 
dz ( Znz — Znz- 1 ) 


( 13 ) 
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is used to approximate the termination plane acoustic pressure gradient. In Eq. (13), p m {z n z, x ) and p m (z n z- U x ) are the 
complex acoustic pressures measured by two lower wall microphones located at z = z, nz - 1 and z = Znz= L, respectively. 
Here nz is the number of axial grid lines in the finite element grid. Just as at the source plane, this acoustic pressure 
gradient profile is also constant across the duct and only the zeroth order horizontal mode (i.e., m = 0) need be 
analyzed. Finally, the measured lower wall acoustic pressure is used as the source- and termination-plane acoustic 
pressure profiles 

p m (0,x) = p m (zi, x i) (14) 

p m (L,x) = p m (z nz , x t) (15) 

where p m (z,\ ,X] ) and p m {z n zi x l) are obtained from the flush-mounted lower wall microphones located at the source 
(i.e., z = 0) and termination (z. = L) planes, respectively. Further, because only plane waves are cut on in the GFIT, 
only the zeroth order horizontal mode (i.e., m = 0) need be analyzed. 


B. Implementation in the CDTR 

Because the CDTR contains higher-order modes that are cut on in the hard wall sections upstream and downstream 
of the test section, the implementation of the PGM to the CDTR is considerably more complicated than the GFIT. 
In the CDTR there are two measurement sections with microphones, one upstream of the liner test section and one 
downstream. 13 The microphone arrays are identical in both sections, each containing 47 flush-mounted microphones. 
The microphones in each array are strategically located around the circumference of the duct at several axial locations. 
The locations are selected to maximize resolution of the modal content of the acoustic field. 13 Data acquired at these 
two microphone arrays are used to obtain the necessary data for the PGM. 

Fet p u (zi, x i,yi) and p d (zi,xj,yj) represent the measured acoustic pressures at microphone I in the upstream (0 < 
zi < Li) and downstream {Li < zi < L) microphone array, respectively. The acoustic pressures at microphone, /, are 
expanded into a series of hard wall duct modes in the hard wall sections upstream and downstream of the liner test 
section as follows 
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Here, NH and NV are the number of horizontal and vertical modes, respectively, in the mode expansion. Generally 
speaking, NH and NV are chosen so that all cut on modes and a couple of cut off modes are included in the series 
expansion. 

Equations (16) and (17) are a set of linear algebraic equations for each microphone location, (z/,x/,y/). If the series 
on the right-hand side of Eqs. (16) and (17) are truncated to the same number of terms as the number of microphones, 
then the two systems of equations can be solved to obtain the mode coefficients of the right and left traveling acoustic 
waves (i.e., A^ mi B^ m ). These known mode coefficients are then used to expand the acoustic pressure field in terms of 
hard wall duct modes: 
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The source and termination plane acoustic pressure profiles ( S m (x) and Q m (x)) needed for the PGM are now determined 
from the mode series 

nv - 1 /nnx\ 

Sm( x )= E [ A mn+ A mn] co ^[-^r ) ( 21 ) 

77=0 
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The measured source and termination plane pressure gradient profiles needed by the PGM are also obtained from these 
mode series 
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V. Results and Discussion 


The primary purpose of the results presented in this section is to validate the PGM for a range of test conditions. 
This is achieved by comparing the educed impedance spectra from the PGM with either known impedance spectra or 
with impedance spectra educed from the WPM. Impedance eduction data were acquired from three primary sources, 
(1) a mode solution that is possible when the liner extends from the source to the duct termination planes, (2) data 
simulated using the finite element solution (i.e., the hard wall sections are included), and (3) data measured in the GFIT 
and the CDTR (i.e., the hard wall sections are included). Three validation exercises are performed. The first validation 
exercise is performed with the CDTR geometry, for which higher-order modes are cut on in the hard wall sections over 
the majority of the frequency range of interest. Results are presented for both a hard wall duct and a conventional liner. 
For the hard wall duct, the known solution is used to simulate the data for impedance eduction. For the conventional 
liner, the normalized acoustic impedance of the liner is selected to correspond to a liner previously tested in the GFIT. 
The eduction data are then simulated from either the mode or finite element solution. Essentially, the first validation 
exercise constitutes a consistency check to demonstrate that the PGM is capable of properly handling the effects of 
higher-order duct modes and flow. Following this consistency check, a second validation exercise is performed. The 
second validation exercise uses measured data from the GFIT (which is capable of acquiring data for both the PGM 
and the WPM) for plane wave test conditions. The third validation exercise uses measured data from the CDTR for 
which higher-order modes are cut on. 


A. The Consistency Check 

The results of the consistency check are restricted to the geometry of the CDTR for which L\ = 0.203 m, Lz = 1 .016 m, 
L = 1.219 m, W = 0.381 m, and H = 0.152 m. As many as three transverse modes may be cut on for the lowest order 
horizontal mode in the frequency range of interest (i.e., 0.300 < / < 2.80 kHz). Results are obtained for a flow off 
(Mq = 0.000) and two flow on conditions (Mo = 0.275, and Mo = 0.500). All results for the consistency check are 
computed at standard atmospheric conditions (cq = 344.28 m/s and po = 1.200 kg/m 3 ). 


1. Hard-Wall Duct Admittance Eductions 


Several examples are now presented to demonstrate the accuracy and robustness of the PGM in hard wall ducts that 

contain planar and higher-order mode sources. Potential problem areas that may require special treatment are also 

highlighted. Because the impedance of a hard wall is effectively infinite, the results are presented in terms of the 

educed normalized admittance. Bo = y- = Ko + ia o- The normalized admittance of the upper wall is set to a zero value 

<=o 

(R// = J = 0.0 + 0.0;), and the normalized conductance, Ko, and susceptance. Go- of the lower wall are educed using 
Qh 

the PGM. The eduction data for these hard wall consistency checks are simulated from the exact mode solution for 
outgoing waves in an infinite hard wall duct: 


p(z,x,; y) 


NH-INV-l 

E E A m, n e iK ™ Z COS 


m = 0 n — 0 



(25) 


Equation (25) is used to extract the source- and termination-plane acoustic pressure and gradient profiles that are 
required for the PGM. Both cut on and cut off modes are allowed to propagate through the flow duct and the normalized 
conductance and susceptance values are educed for the following four sources 

1. a plane-wave source (i.e., Aq 0 = 20 Pa) 
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2. the (01) mode (i.e., Af t = 20 Pa ) 

3. the (10) mode (i.e., A^ () = 20 Pa ) 

4. the (11) mode (i.e., 4 jj = 20 Pa ). 

Prior to performing the impedance eductions, a grid refinement study was conducted using each of the four sound 
sources to determine the coarsest finite element grid for accurate resolution of the acoustic waves in the hard wall duct. 
The finite element grid used to perform the impedance eductions was four times denser than the coarsest grid needed 
for accurate resolution of the acoustic wave in the hard wall duct. 

Figure 3a shows the educed normalized conductance and susceptance for the plane wave source (i.e., the (00) 
hard wall mode) for the flow off ( Mq = 0.000) and the two flow on conditions (Mq = 0.275, and Mq = 0.500). The 
conductance and susceptance were educed every 0.100 kHz for frequencies ranging from 0.300 to 2.800 kHz. As can 
be seen from the figure, the normalized conductance and susceptance values that were educed with the PGM are in 
excellent agreement with their expected value of zero. There are 144 test conditions (i.e., three Mach numbers and 

48 frequencies) and approximately 10% of these test conditions had to be restarted because the SDFP did not give 

a normal return using the initial starting value of Ko = Oo = 0.500. In each of these cases the initial starting value 
had to be reduced below 0.500. Figure 3b shows the educed normalized conductance and susceptance for a higher 
order mode (i.e. the (01) hard wall mode). The cut on frequency for this mode is approximately 1.200 kHz. At 
frequencies below the cut on frequency, the sound wave is rapidly attenuated, whereas at frequencies above the cut 
on frequency the sound wave propagates without attenuation. Several observations may be made from Fig. 3b. First, 
the conductance and susceptance at all frequencies below the cut on frequency are accurately educed. Second, the 
conductance at frequencies above the cut on frequency is also accurately educed. Third, there are some errors in the 
educed susceptance at a few of the frequencies just above the cut on frequency. This trend has also been observed 
in studies with the WPM 7 and is due to the fact that above the cut on frequency the sound wave is not attenuated. 
This leads to an objective function that is extremely flat at these frequencies, causing the educed susceptance to be 
highly dependent on the starting values. Results for the (01) and (11) hard wall mode are not presented for the sake 
of brevity. However, the educed conductance and susceptance values for each of these higher-order modes produced 
trends similar to those observed in Fig. 3b. 


2. Soft -Wall Duct Impedance Eductions without Hard Wall Sections 

The results in this subsection test the capability of the PGM to accurately educe the soft-wall impedance. The PGM 
impedance eduction data were obtained from a mode solution that is possible when the upper and lower wall liners are 
uniform and cover the full extent of the duct (i.e., extend from z = 0, to z = L). In this situation the quasi-3D mode, 
p m {z,x), has the following closed form solution 

Pm{z,x) =X mn (x)e lK ™ z (26) 


where n is used to order the modes according to their attenuation rates. The vertical modes, X mn (x), are obtained using 
separation of variables 


j cos (X„x) + i !i a j, nn sin (X„x) , for X n ± 0 

Xmn(x) — ^ n ^t) 

1 , for X„ = 0, or a„„, = 0 


where X n are complex zeroes of the transcendental function 

' 1 

C h ) UoU 


E(k,H,MoXoXH,Xn) = ikHa mn X n H I - ( \^d} mn + X\H 2 ) tan (X„H) 
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Note that for a given value of X„ the axial propagation constant, K m „, that are the complex zeroes of the trancendental 
function, E(k,H,Mo,C,oXH,X n ), has two roots. It is easily shown that right moving waves (i.e., K+ n ) are identified 
as those modes for which the axial propagation constant has a zero or negative imaginary part and positive real part, 
respectively. The second root (i.e., K mn ) corresponds to left moving waves. In a hard wall duct where, A„ = , 

Eq. (30) reduces to the expression given in Eq. (18). 

In this paper we are interested in the aft-duct configuration (the sound and the flow are propagating in the same 
direction) so only Kf n is of interest. The method of analysis used here calculates the mode eigenvalues, A„, that are 
the complex zeroes of the transcendental function, E(k,H,Mo,C,o,C,H,X n ), and the corresponding axial propagation 
constants, K^ n . The method used to obtain these roots is an analytic continuation method similar to that used by 
Snider, 20 where the zero-flow modes are calculated using a matrix method and the flow modes are calculated from the 
zero-flow modes by gradually increasing the mean flow Mach number from its zero value. 

Given the quasi-3D vertical mode, X mn (x) , and axial propagation constant, K nun the closed form solution in Eq. (26) 
is used to extract the PGM impedance eduction data in the aft-duct configuration for a specified horizontal mode (i.e., 
m) and a given vertical mode, n 


Sm (x ) — X mn (x) , Q lr i (A ) — X n 


n (x)e‘ K " mL , 


dp(0,x) 

dz 


— iK+ n X mn (x ) , 


dp(L,x) 

dz 


= iK+X„ 


,{x)e lK " mL 


(31) 


The PGM impedance eductions are performed for the first two vertical modes (n = 1 and n = 2) where the modes are 
arranged according to their attenuation rate. The test liner is a single-layer, conventional perforate-over-honeycomb 
liner. This liner was chosen because it is representative of liners commonly used for the reduction of fan noise in 
current aircraft engine nacelles. The facesheet for this conventional liner has 0.991-mm-diameter holes, an open area 
ratio of 0.087, and a thickness of 0.635 mm. The depth of the cavity is 38.18 mm; this depth was chosen so that the 
resonant frequency of the liner is in the middle of the frequency range of interest. The impedance of this liner was 
determined from earlier tests in the GFIT and these impedance values were used in the transcendental equation to 
extract the vertical modes and impedance eduction data for the PGM. 

Figures 4a and 4b compare the normalized resistance and reactance values for the (01), (1 1), (02), and (12) modes 
educed by the PGM with that measured in the GFIT for the flow off condition. Results are presented in 0.200 kHz 
increments. For each sound source, the PGM educed resistance and reactance values are not distinguishable from those 
obtained using the WPM in the GFIT. Note that when the flow is off, the conventional liner is a low resistance liner 
whose resistance is nearly independent of frequency (as expected), and the normalized reactance follows a — cot (kd) 
(where d is the depth of the cavity) behavior with the resonant frequency (i.e., the frequency at which the reactance 
approaches zero) at 1.500 kHz (as expected). 

The normalized resistance and reactance spectra educed with the WPM and the PGM for a flow on condition (Mach 

0.275) and several higher order mode sources (the (01), (11), (02), and (12) modes) are depicted in Figs. 5a and 5b. 
Educed impedances from both methods (i.e., the WPM and PGM) are in good agreement with each other. The PGM 
educes the same impedance for each higher order mode source as expected. Both eduction methods show an increase 
in the normalized resistance as the flow Mach number increases from zero (as expected), but the reactance at the low 
end of each spectra varies from the — cot (kd) behavior that was observed for the flow off condition. The differences 
at the low frequency end of the spectrum were expected and result from the fact that the frequency is approaching an 
anti-resonant value. This result corresponds to previously reported findings in which the reduced attenuation due to 
an impedance near anti-resonance causes the “eye” of the objective function contour to be quite diffused (i.e., a large 
range of impedances produces similar attenuations). 


3. Soft-Wall Duct Impedance Eductions with Hard Wall Sections Upstream and Downstream 

Results presented thus far test the capability of the PGM to accurately educe the conventional liner impedance when 
a closed-form solution is available to extract the impedance eduction data. However, in real facilities the flow duct 
contains hard wall sections upstream and downstream of the liner test section so that a closed-form solution is not 
available. Therefore to test the PGM in the presence of more realistic test conditions, simulations are performed in 
the CDTR geometry with the hard wall sections included. The axial pressure gradient that is needed to educe the 
impedance is obtained as follows: 

1 . The source- and termination-plane acoustic pressure profiles for the CDTR are set to the same profiles that were 
used in the hard wall example problem. 

2. The finite element method is then used (with the conventional liner impedance that was educed in the GFIT) to 
simulate the axial acoustic pressure gradients needed for the PGM. 
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PGM impedance eductions with hard wall sections upstream and downstream of the conventional liner were computed 
with the flow off (Mach 0.000) and flow on (Mach 0.275). For each sound source (i.e., the (00), (01), (10), and (11)), 
the normalized resistance and reactance values educed with the PGM matched those educed in the GFIT to two 
decimal digits of precision except at 0.4 kHz where the attenuation is low. Graphical results for these test cases are not 
presented in this paper for the sake of brevity. 

B. Validation Using Measured Data from the GFIT 

PGM results are now compared to WPM results for cases in which measured data from the GFIT is used. Results 
are presented for both a hard wall insert and the conventional liner. Data were acquired in the GFIT with the flow off 
(Mach 0.000) and for two flow on conditions (Mach 0.300 and 0.500). The Mach number for the flow-on condition is 
the centerline value. The GFIT geometry is rectangular with H =0.060 m, W = 0.050 m, L\ = 0.200 m, L 2 = 0.610 m, 
and L = 1.220 m (see Fig. 1). For the WPM, measured acoustic pressure data were obtained at 53 unevenly spaced 
points along the bottom wall. For the PGM, a measurement of the acoustic pressure is taken only at the source and exit 
planes and at 0.030 m away from these planes (i.e., only four acoustic pressure measurements are required). Because 
only plane waves are cut on at the frequencies of interest, all acoustic pressure measurements are obtained at the lower 
hard wall, and these are assumed to be constant in the vertical direction of the duct (i.e., from the lower wall to the 
upper wall). The hard wall insert is a 12.700-mm-thick stainless-steel plate and is chosen because it provides a liner 
for which the normalized admittance, [)//, is known (i.e., [)// = 0.000 + 0.000/ for a hard wall insert). 

1. Hard Wall Duct Admittance Eductions 

Figure 6 compares the educed normalized conductance and susceptance spectra for the hard wall insert for the flow-off 
condition. Both the WPM and the PGM are in excellent agreement with each other and with the known conductance 
and susceptance for a hard wall insert. Figure 7 shows the comparisons for the two flow-on conditions. The agreement 
between the normalized conductance spectra educed with the WPM and the PGM is quite good for Mo = 0.300 but 
degrades at Mo = 0.500. The normalized conductance is observed to be fairly close to the expected value of zero for 
all conditions, but the normalized susceptance deviates from the expected value of zero for both the WPM and PGM 
for the flow on conditions. This scatter in the susceptance spectra about the expected value of zero increases with 
the mean flow Mach number. This result is not surprising when one considers that the refractive effects of the mean 
boundary layer (which are neglected in both the WPM and the PGM) are more important at higher mean flow Mach 
numbers. 

2. Conventional Liner Impedance Eductions 

Figure 8 compares the educed normalized resistance and reactance spectra for the conventional liner for the flow-off 
condition (i.e.. Mo = 0.000). Both the resistance and reactance spectra from the WPM and the PGM method are in 
excellent agreement with each other. The normalized resistance and reactance spectra educed with the WPM and the 
PGM for both of the flow-on conditions are depicted in Fig. 9. Educed impedances from both methods (i.e., the WPM 
and PGM) are in good agreement with each other, especially at the lower Mach numbers. 

C. Results Using CDTR Measurements 

The PGM results for the CDTR were also computed. These were performed as follows: 

1. Incident sound waves of 130 dB (with horizontal mode order, mo) were generated upstream of the liner test 
section 

2. The sound field at each microphone array in the CDTR was simulated using 20 modes propagating upstream 
and 20 modes propagating downstream (i.e., NV = 5 and NH = 4 in each hard wall section). With this number 
of modes the microphone acoustic pressure residual (i.e., the difference between the measured and simulated 
microphone acoustic pressure) was no larger than 1 0 1 0 

3. Equations (21)-(24) are used to acquire impedance eduction data using m = mo. The quasi-3D equation (Eq. (2)) 
with m = niQ is used in conjunction with the SDFP algorithm to educe the impedance of the liner 

Figure 10 shows the educed normalized admittance for the plane wave (i.e., the (00) mode) and the (01) mode for 
the flow off condition. With few exceptions, the educed conductance is observed to be close to the expected value 
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of zero for each of these source conditions. The conductance values educed with the (00) mode deviate from zero 
at 1.5 kHz, while those educed with the (01) mode show scatter at 1.1 and 1.5 kHz. On the other hand, the educed 
susceptance exhibits significantly more scatter than observed for the conductance, for both the (00) and (01) modes. 
Thus, although the predicted attenuation remains low, the current implementation of the PGM is clearly insufficient for 
this configuration. All indications are that this problem is caused by the mode expansion at the microphone locations. 
One solution may be to use the measured acoustic pressure data directly, avoiding the use of the modal expansion. 
However, this approach has the disadvantage that the solution to the BVP would have to be obtained with a 3-D code 
because higher order modes are cut on in the CDTR. Potential solutions to this issue, along with other improvements 
to the PGM, are currently being explored. 


VI. Concluding Remarks 

The pressure gradient method (PGM) of impedance eduction has advantages over the wall pressure method 
(WPM). The primary one is that it does not require the installation of microphones in the liner test section to educe the 
impedance. This paper has tested the PGM over a range of test conditions using eduction data that were (1) obtained 
from a closed-form solution, (2) synthesized from a finite element method, and (3) measured in the NASA Langley 
Grazing Flow Impedance Tube (GFIT) and Curved Duct Test Rig (CDTR). Based upon results presented in this paper 
the following primary conclusions are drawn: 

1 . When used with the closed-form and synthesized eduction data, the PGM educed impedance converges almost 
exactly to the known impedance of a hard wall and to that of a conventional test liner 

2. When used with measured data from the Langley Grazing Flow Impedance Tube (GFIT), the PGM produces 
results that are consistent with that of the Wall Pressure Method (WPM) 

3. Validation of the method in the Curve Duct Test Rig (CDTR) is continuing and may require the direct use of 
measured data without mode decomposition. Multiple options are being considered to support this approach 
with quasi-3D analysis. However, it appears likely that a fully 3-D code will be required to achieve robust 
impedance eduction in this aeroacoustic environment. 
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Figure 1. Schematic of the GFIT with microphones embedded in the liner test section. 
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Figure 2. Schematic of the CDTR flow duct with microphones removed from liner test section. 
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Figure 4. PGM educed normalized resistance and reactance for conventional liner for several higher order mode sources with flow off. 
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Figure 5. PGM educed normalized resistance and reactance for conventional liner for several higher order mode sources at Mach 0.275. 
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Figure 6. Educed normalized conductance and susceptance of the hard wall insert for the flow-off condition. 
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Figure 7. Educed normalized conductance and susceptance of the hard wall insert for the flow-on conditions. 
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Figure 8. Educed normalized resistance and reactance of the conventional liner for the flow-off condition. 
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Figure 9. Educed normalized resistance and reactance of the conventional liner for the flow-on conditions. 


6 r 


A Educed normalized conductance (00) mode 
V Educed normalized susceptance (00) mode 
> Educed normalized conductance (01) mode 
- <] Educe normalized susceptance (01) mode 


ID 

o 

C 

(0 

E 

T3 

CO 

T3 

CD 

N 

75 

E 

o 


2 - 


-2 - 




0.5 




1.5 

Frequency, kHz 


Figure 10. Educed normalized conductance and susceptance using the PGM for a hard wall insert in the CDTR with flow off . 
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